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I. INTRODUCTION 



At a temperature much lower than the particle mass, it is physically reasonable to expect 
that the leading order thermal correction to the physical self-energy be given by the forward 
scattering amplitude integrated over the thermal phase-space of the lighter particle. For 
instance, in the literature, the thermal self-energy of a particle is often written as fl|- 



d 3 k 

Zr(l>) = I (27r) 3 2 £ ^ste) T(k + p^k + p) , (1) 



where -E^ght = yk 2 + m 2 ight is the energy of the lighter particle with mass m light , n(Ei ight ) 
is a Bose or Fermi distribution function, and T(k + p — ► k + p) is the forward scattering 
amplitude related to the usual / by 

T(k+p-^k+p) = 8ir v Gf{k,p) . (2) 

The minus sign on the right hand side of Eq. ([!]) stems from the fact that the scattering 
amplitude is defined to be —i times an iV-point correlation function whereas the self-energy 
is defined to be i times a 2-point correlation function. 

A merit of this method is that one can simply take the scattering amplitudes from 
experiment. In reactions involving strong couplings, this may be the only reliable way of 
calculating the thermal correction to the real part of the self-energy [f|]. As the density 
and temperature become higher, Eq. ([]]) needs higher order corrections. In this paper, we 
start with the imaginary-time formulation of relativistic finite temperature field theory and 
obtain an exact expression for the thermal self-energy. 

There are a variety of finite temperature field theory methods used in literature. The real- 
time method was first investigated by Schwinger and Keldysh and later in a slightly different 
and modern form by many people |5|-|9|]. The imaginary-time method was introduced by 
Matsubara ]TD] and through analytic continuation can be related to the real-time method 



TT| . Of course any equilibrium Green function can be calculated using any of the methods, 
given the analytic relations between them. 

There are, however, situations where one particular method is more economical. For 
instance, if one would like to study time-ordered real-time correlation functions, the real- 
time method is very much the natural choice. If one would like to study a response function 
(a retarded function), the most economical way would be to compute the imaginary-time 
correlation function with Matsubara frequencies and then analytically continue the result. 

In this paper we are interested in the retarded correlation functions. We would like to 
use the imaginary-time method. However, since the imaginary time is a fictitious parame- 
ter introduced to deal with the temperature, it is not easy to gain physical insight by just 
looking at the end result of the imaginary time calculation. A way to remedy this difficulty 



was developed by one of the present authors Jl2| based on earlier works (T^-|T5[. There, a 



diagrammatic method to calculate the spectral density of two-point functions in a scalar 
theory was presented, starting from the imaginary-time formulation of finite temperature 
field theory. Subsequently it was used to calculate the leading order hydrodynamic coeffi- 
cients in a scalar theory [l(J. In this paper the above work is extended to the calculation of 
iV-point retarded functions. 

The analysis in this paper is based on the fact that applying the analytic continuation 

iu -> k° + ie (3) 

to each of the independent external Matsubara frequencies produces the (Fourier trans- 
formed) retarded correlation function |5,17|. When the chemical potentials all vanish, this 



operation produces the physical retarded function. However, with non-zero chemical poten- 
tials one must be more careful. 

In the imaginary-time formalism, finite density is dealt with by using an effective Hamil- 
tonian 

K = H-^2^aQa, (4) 

a 

where fi a is the chemical potential associated with each conserved charge Q a - Hence, the 
retarded function produced by using Eq. (|3|) corresponds to the retarded correlation function 
of <pK{t) = e tpe~ . But the physical retarded function is a correlation function of 
ifH(t) = e tHt 0e~ lHt , not Fortunately, there is a simple relation between the two. In 

most cases the fields of interest are charge eigenstates. Then 

K {t) = e^W*) , (5) 

where 

V<p = ^ a ( 6 ) 

a 

is the total chemical potential associated with the field <p carrying conserved charges q a . In 
momentum space, the two retarded functions are related by 



OW) = 0fc-M), (7) 

where we used a shorthand notation, k — /i v = (A; — fi^, k). As is shown in Sec. [II B| , this 
shift also has the effect of separating the purely dynamic part of the theory and the purely 
statistical part. 

This paper is organized as follows. In Sec. |H] a brief derivation is given for the diagram- 



matic rules to compute an arbitrary retarded function. Section |TJ presents a derivation of 
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the self-energy formula using the results from Sec. 0. As an example, the self-energy of an 
electron at the one-loop level is worked out in Sec. [IV]. Finally, we conclude in Sec. |V[ Ap- 
pendices [A] and [B] list some relevant facts about the finite temperature propagators used in 
this paper. Appendix |C| deals with derivative couplings. Appendix [D] amplifies the discussion 
of the thermodynamic potential in Sec. [Hi] . 



II. iV-POINT RETARDED CORRELATION FUNCTIONS 



In this section, as a prelude to the self-energy calculation, diagrammatic rules to calcu- 
late an arbitrary A-point retarded function are derived starting from imaginary-time finite 
temperature field theory. The derivation here closely follows Ref. []T^[. Details omitted for 
the sake of brevity and readability can be found in that reference and references therein. 
The rules derived here coincide with the rules obtained by Evans || who used the real-time 
method and the analytic relations between A-point functions. For simplicity, only theories 
involving no derivative couplings are considered here. However, the argument given here 
can be straightforwardly generalized to theories with derivative couplings, such as chiral 
perturbation theory, as indicated in Appendix £]. 



A. Euclidean Correlation Functions 



It is known that performing Matsubara frequency sums in an Euclidean Feynman diagram 



produces terms corresponding to old-fashioned time-ordered perturbation theory |]14 |. The 
aim here is therefore to provide conventions and definitions along with a sketch of a proof, 
but not the details. A detailed derivation for a relativistic case can be found in Ref. [T2| 



and a non-relativistic version can be found in Ref. [13| . 

The standard Feynman rules use the momentum space Feynman propagator. For our 
purposes, it is more convenient to use mixed propagators which are functions of time and 
spatial momentum. As shown in Appendix A, the mixed propagator for both bosons and 
fermions can be represented by 

G c (r,k) = J ^ A C H ( e —p+(*)0(r) +e«Tp-(k)9(-T)) . (8) 

Here u = /c°, and the spectral densities are labeled pf(k) or p^(u,k), as convenient. The 
label ( distinguishes boson or fermion, {B, F}. The statistical factor N^(u) is 

%) = i + (-i)SH. (9) 

where (— 1) B = 1 and (— 1) F = —1 and 
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"cM = ^ - (_i)c • ( 10 ) 

Explicit forms for the spectral densities are not of interest in this section. What is important 
for us is that they satisfy 

p+(u,,k) = -^(-u;,k) (11) 

due to the periodicity or antiperiodicity of G^ (r, k) in r. A few relevant facts about the 
spectral densities are given in Appendices [A] and [TJ. 

The Feynman rules in this mixed space of imaginary-time and momentum are al- 
most identical to the standard momentum space Feynman rules. The differences are: An 
imaginary-time Tj labels each vertex including the external ones. Each line connecting two 
vertices labeled by the times r a and r& represents a propagator G^(r a — r b , |k|) if the mo- 
mentum flows from t& to r a . If the momentum flows from r a to r&, the line represents 
G^(Tb—r a , |k|). The direction of the momentum should follow the direction of charge flow. 
Instead of the sum over all loop frequencies, there are integrations over all Tj's from to (3. 
At each vertex where an external field operator extracts external frequency iuj an additional 
factor of exp (iVjTj) is present. Hence, the contribution of a connected Feynman diagram 
T with a total of V+l vertices and N+l external operator insertions has the following 
schematic form: 

C ( P +1 ({qi,iui}) 

(12) 

Here the independent spatial loop momenta are denoted by k^, and the quantity Ay+\ 
includes all the factors from interaction vertices such as coupling constants, 7 matrices, and 
symmetry factors. We suppress Lorentz indices, if there are any, since they are not relevant 
for the moment. A given line in the Feynman diagram running between vertices acting at 
times r a and r& is denoted by a. 

We know that as the temperature goes to zero the result of the time integrations must 
be that of old-fashioned time-ordered perturbation theory. That is, 

KmC i p +1 ({q 1 ,ivi})= £ hm^ r ;!({q,,^}), (13) 

where a represents a given permutation of the vertices and T a represents a time-ordered 
diagram for that permutation given by 
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limCi r ;i({q ; ,^}) 



d*k L TT ( f e{Ua) pt{ka \ A * t l JJ (AJ-uf)" 1 . (14) 




/ intervals 
V>j>l 

Here we omit the overall frequency-conserving Kronecker-5. The sign s a = {+, — } depends 
on both the time ordering and the direction of the charge flow. If the charge flows in the 
same direction as the time a + sign is assigned, otherwise s a = —. If there are no conserved 
charges the distinction is without significance. The prefactor Ay+\ is the same as before. 
The time ordering also determines the frequency denominator (Aj — ivj)- Given a time 
ordering a, AJ is given by the sum of frequencies of all lines crossing the j-th interval IJ 

AJ = 5>« • (15) 
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J 



Similarly, uj is the net external frequency flowing out of the diagram above the interval IJ 



»■= E ( 16 ) 

V>l>j 



In this limit the external frequencies v\ are continuous. 

Comparing Eq. (|14|) and Eq. (|l"2D one sees that one should keep all factors of N^(u c 



since 



lim N c (u a ) = 9{u a ) . (17) 

p— >oo 



Then for Eq. ([14]) to be true, there cannot be any additional factors of exp {±/5 J2i w «} m 
the finite temperature result since they will make the result diverge due to the fact that 
frequencies can be both positive and negative. The structure of integrand in Eq. (|T2|) also 
dictates that vanishing contributions such as exp{— f3\ Y^i^iW or l//5 n cannot occur. Then 
the only feasible non-zero temperature result is Eq. ( |j"4j ) with 9{ui a ) replaced by N^ a {uj a ). 
That is, 



cSU{qiM})= E C^\({q h m}) , (18) 

where 



r CT cr 



CK({qj,^}) 

=/n n (/ ^jvc-w/sew) <s . n (aj-<) _i , <i<» 



LgT \ / o6T \ / intervals 

V>j>l 



(2 

again omitting the overall frequency-conserving Kronecker-5. 
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B. Retarded Correlation Functions 



The (iV+l)-point retarded function is denned by |L9 

R(x] {xi}) 

x 

II ^(^ J -^+i)([---[[^( ;r )>^ 1 (a; (71 )],^ (T2 (x (T2 )],...,^ (TJV (x ffiV )]), (20) 

a j = l 

where t = x° is the fixed largest time, % runs from 1 to N, and Y^a indicates the sum over 
all possible permutation of {£i,i2, • • • , tw}- The operator if can be any field in the theory. 
When two fermionic quantities are involved, the commutator becomes an anti-commutator. 
The angular bracket represents the thermal average. From here on, we always denote by v 
the vertex at fixed time t. 

Several authors |5|jr?J] have shown that the following analytic continuation of the (N+l)- 
point imaginary-time correlation function leads to the retarded function: 

ivi — > qf + ie for I ^ v 

ivv -^q° v + iNe , (21) 

where at vertices other than the fixed time vertex v the frequencies are injected, and at v 
the frequency is extracted so that the frequency is conserved including the small imaginary 
part: 

q° v + iNe = + ie) . (22) 

The contribution of a Feynman diagram T to the retarded function R y N ^\({qi}) (subscript v 
signifies that v is the fixed vertex) is then: 



r CT cr J Ler^' 1 ) aer V z/l / 

x II (AJ-^ + ie)" 1 n (Af-Vf-ie)' 1 . (23) 

intervals intervals 
V>j>v+1 v>j'>l 

The AJ were defined in Eq. (|T5D and the quantities qj refer to the net external frequency 
flowing out of the diagram above the given interval 

Qj= E (24) 

V>v t >j 

Here we chronologically label the vertices by {to, • • • ,ty} and the interval between tj and 
tj_i is the j-th interval. Also, a product without a factor is defined to be 1. Note the change 
in the sign of the e term when crossing the v-th interval. The sign depends on whether 
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the interval is above or below the vertex v which extracts iNe. Intervals where no external 
frequencies flow in or out have no ie to start with, but it is not hard to make the internal 
frequencies slightly off the real axis to make a consistent assignment of the imaginary parts 

0- 

We would like to express Eq. fl23|) in terms of uncut and cut propagators in analogous 
fashion to the Cutkosky approach at zero temperature. To that end, notice that if all the 
+ie's in the frequency denominators were — ie and all the N^ a (cj a ) were 9(u a ), then the 
expression would be exactly that of old-fashioned real-time perturbation theory. When the 
time orderings are summed, this expression just reproduces the Feynman diagram expression. 
Hence when all the small imaginary parts are — ie, all we need to do to sum the time-ordered 
terms is to change the zero temperature propagator 

G» W = / £ ( + M.) (25 ) 
s J 2m \uj — k v — ie to + — ie) 

to the finite temperature one 
see Appendix B. Then 

- E r / n ^ n (/ ^w^ flI ik Q |)) A»ti 




(2vr) 4 



^vSIIW- ( 27 ) 



Lev ) «er 

The factor — i(i) v+1 comes from having V integrals with — ie. The phase i v+1 can be 
absorbed into Ay+\ by adding a factor of % to the coupling constant at each vertex as in the 
usual diagram rules, e.g. see Peskin and Schroeder |20) whose conventions we follow. 

We would like to apply this argument to each of the ±ie parts of the product in 
Rn+i({1i})- in order to do so, the two parts must be separated from each other so that the 
lines common to both parts can be regarded as external ones. This is achieved by using the 
identity 

o= n {uj-iey 1 - n (V^r 1 -* e n fa+uy 1 n 



, 3 -ie 



j=v+l j=v+l £=«+! fc=(+l j=v+l 



(28) 



Here all U/s are real, a sum without a summand is and a product without a factor is 1. 
This identity follows directly from 2iri5(x) = (x — ie)^ 1 — (x + ie)^ 1 . We refer to the lines 
corresponding to the interval in 5{U{) as the cut lines. 



7 



Using the above identity and performing the relevant resummation of time ordered terms 
on both sides of the cut, we obtain 

R { nI\({h}) = D ( P +1 ({ qi }) + £ J$;\({«}) , (29) 

where 

RpMiri) = -i{-i) v+ ar- 1 n kxi n pz%)) 

LeT V / cut lines 

c 

x n G c a (^)* n gum, (3°) 

«er + «'er_ 

and the sum is over the cut diagrams ¥ v . The uncut propagator corresponds to the lines 
in the side T_ which contains v, and G£ corresponds to the lines in the other side r + . For 
the cut lines the direction of time flow is from T_ to T + . Here V± denotes number of vertices 
in T±. Again, the factors (— i) v+ (i) v - can be absorbed into A v ^\ by including a factor of i 
in the contribution of each vertex in T_ and a factor of —i in the contribution of each vertex 
in r + . 

So far, nothing has been dependent on the presence of chemical potentials. When 
chemical potentials are present, Eq. ( |29| ) corresponds to the retarded function of fields 
<fK(t) = e lKt 0e~ tKt with K = H — J2 a paQa- As explained in Sec. |, there is a simple 
relation between (p K and <pn, <pi((t) = e lfJ,ipt (pH{t)- Hence, to convert a retarded function of 
the (fix to the corresponding retarded function of the 0h, all one has to do is to shift each 
external frequency from q to q—^ip if it corresponds to the particle of the species and to q+^Lp 
if it corresponds to the anti-particle. These shifts have important simplifying consequences. 
Due to charge conservation, not only energy-momentum, but also chemical potential has 
to be conserved at each vertex. Hence, shifting the external frequency q° — > q° — fi^ is 
equivalent to shifting all internal frequencies 1° — > 1° — fi a according to the species a. 

The significance of these internal frequency shifts is that they remove the chemical po- 
tential dependence from the spectral densities. For instance, when the chemical potential 
is non-zero, the spectral density of the bosonic propagator appearing in Eq. Q23D is (see 
Appendix [S]) 

pUk°) = sign(fc° + ,x)27r 5 ((k° + rf- El) , (31) 

where Ek = -y/k 2 + m 2 . By shifting k° — > k° — [i, p^(k° — p) becomes independent of /a, 
while the same change will make the statistical factor become n B (u — /i).[| 



1 These internal shifts also have a simplifying effect on derivative couplings. The derivative dt 
present at an interaction vertex yields 1° + fj, which becomes simply 1° after the shift of internal 
frequencies. 
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For small chemical potentials such that > even more simplification is possible 
because in that case, 9{E ± /i) — 6(E) and 6(—E ± /i) = 9(—E). Then 

iV c (A; - //) = 9(k°) + (-l) c sign(£; )n c (|fc - /i|) , (32) 

when k° = ±-E&. Consequently, the chemical potential only appears in the statistical factors. 
In this way we achieve the separation of the statistical part and the purely dynamic part of 
the particle propagation. (For a slightly different approach, see PTJ). 

We can now state diagrammatic rules for calculating an arbitrary iV-point retarded 
function with small chemical potentials (|/^| < m^), see Appendix [B]for further details. For 
convenience, we denote the thermal phase space factor for scalar bosons by 

r B (fc) =n B {\k° ~ fi\)27T5{k 2 -m 2 ) (33) 

and for spin-| fermions 

T F (k) = -n F {\k° - / u|)(P7 M + m)2vT(5(A; 2 - m 2 ) . (34) 
Define the cut propagators for scalar particles 

A±(k) = 6(±k°) 2tt 5(k 2 - m 2 ) + r B (Jfe) , (35) 

and for fermions, 

A^(k) = 6{±k°) 2tt (P 7m + m) 5{k 2 - m 2 ) + T F {k) . (36) 
For scalar particles, the uncut propagators are given by 

G B (k) = k2 _l 2 + ie + r B (*) and GUk) = k2 _^_ ie + r B (*) (37) 
and for spin-^ fermions, they are given by 



G F ykj — —z p +r F (fc) and G F (k) = —— <=- — + T F (k) . (38) 

k — m z + te k A — m z — le 

For gauge bosons, set m = in each of the scalar particle propagator, and multiply by 
{—guv) (Feynman gauge). Note that all of the cut and uncut propagators can be written as a 
zero-temperature, zero-density part plus a common thermal phase space factor. Not surpris- 
ingly, the zero temperature parts of the propagators all coincide with the zero-temperature 



Cutkosky rule propagators [22 



We call the region where the fixed vertex v sits the unshaded region and the other half 
the shaded region. The rules to calculate the retarded functions of (pu are: 
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1. Draw all topologically distinct cut diagrams including totally uncut ones, keeping the 
largest time always on the unshaded side. Disconnected pieces produced by the cutting 
are allowed. 

2. Assign momenta to the lines according to the flow of the conserved charges, if there 
are any. 

3. Use the usual Feynman rules for the unshaded side assigning G^(k) to the uncut lines. 

4. Use the complex conjugate Feynman rules for the shaded side assigning G^(k) to the 
uncut lines. The 7 matrices are not to be complex conjugated. 

5. If the momentum k of a cut line crosses from the unshaded region to the shaded region, 
assign A£(k). 

6. If the momentum A; of a cut line crosses from the shaded region to the unshaded region, 
assign A^(k). 

7. Divide by the symmetry factor if applicable. There is an overall factor of —i. 

Equation fl29"|) and the above diagrammatic rules are the main results of this section. As 
shown in Appendix |C], this same set of rules also apply when there are derivative couplings. 

As an example of applying the rules, consider a scalar theory with a —g<f) 3 interaction in 
the Lagrangian. The diagram shown in Fig. p] yields the expression 

C^\k) = {-i)(-ig)\igf j tJ^^ 7^ A B(0AMA B (g) A^l+p-q) 

x G B (l+p-q-k) G B (l-k) G B (p-q) G* B (1+P) • (39) 

Note that at T = 0, the contribution of this diagram is zero because the cut part then 
represents a process where four physical particles annihilate each other into the vacuum. 
At non-zero temperature, this diagram gives a non-zero contribution because it can also 
represent scattering between physical particles in the thermal medium. 

III. MULTIPLE SCATTERING EXPANSION OF THE SELF-ENERGY 

In this section, we present the multiple scattering expansion for the thermal contribution 
to the self-energy. The retarded self-energy is, of course, a retarded 2-point function. From 
the previous section (see also Appendix |TJ) we know that as long as m a > \/J, a \, all propagators 
satisfy 

D c (k) = D°(k) + r (; (k) (40) 
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FIG. 1. A typical finite temperature cut diagram in a scalar gcfft theory. Each cut line con- 
tributes Ag or Ag according to its orientation. In the unshaded (shaded) region, an uncut prop- 
agator contributes Gb (G b ), and an interaction vertex contributes a factor of — ig (ig). 

where r^(fc) is the thermal phase space factor common to all four cut and uncut propagators, 
and D®(k) is the zero temperature, zero density propagator. 

Using Eq. (|40| ) we can now expand the expression (^) for the self-energy in the number of 
T^'s. The coefficients in this expansion involve only the zero temperature propagators. Hence, 
the coefficient function is either a zero-temperature Feynman diagram or a zero-temperature 
Cutkosky diagram. External momenta for this coefficient function are provided by the self- 
energy momentum, k, and the thermal particle momenta from the I^'s. 

We would like to relate the coefficients of the expansion to the physical zero-temperature 
scattering amplitudes. In order to do so, we need to ensure that in the expansion we consider: 

1. Symmetry factors are all correctly accounted for. 

2. Disconnected parts cancel. 

3. Self-energy insertions do not cause divergences. 

4. Additional polarization factors needed for the external thermal particles are all cor- 
rectly provided. 

5. When fermions are involved, the overall sign for an individual Feynman diagram is 
correctly produced. 

We present the result of the expansion first and deal with the above points later in Sub- 
sees. pTAh^TTC . 

We emphasize that the formal expansion in the number of T^'s can be always made. To 
be useful, however, the series must be truncated. When the temperature and the density 
are low so that T <C \m a ± fi a \ for all particle species a, one may order the expansion with 
respect to the relative strengths of e~^ ma ~^ a ^ T and e~^ ma+IJ ' a ^ T . The first few terms of such 
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(a) (b) (c) 

FIG. 2. Diagrams that contribute to the self-energy. The upper row is illustrated in more detail 
in the lower row. 

a virial expansion will be a good approximation as long as the densities and the temperature 
stay low. As the densities and the temperature grow so that the minimum of \m a ± fi a \ is 
no longer small compared to the temperature, we lose control of the approximation because, 
potentially, an infinite number of terms in the expansion become important. 

The diagrams that provide the thermal contributions to the retarded self-energy, Sj; et , 
can be separated into three groups, as depicted in Fig. Q Group (a) contains no cut, group 
(b) contains cuts that do not separate the two external vertices and group (c) contains cuts 
that do separate the two external vertices. 

Consider first the diagrams depicted in Fig. 0(a). The diagrams in this group are solely 
made of the uncut propagator G^(l) = G®(1) + where G®(1) is the zero temperature 

Feynman propagator and r^(7) is the thermal phase space factor. When the expansion in the 
number of is made, the coefficients are obviously the zero temperature Feynman diagrams 
whose external momenta are supplied by the self-energy momentum and the thermal particle 
momenta from the T^'s. That is, the coefficients correspond to matrix elements of the 
scattering operator T, with the 5-matrix given by S = 1 + iT . Hence, the expansion of the 
first group of the diagrams can be written as 

i . n-l 

4 a) (*) = - e «— / n^nu^immun 

n>2,a b W} J i=l 

+ (disconnected diagrams) , (41) 
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where 

d 3 l- 

^•^(^k^-^- (42) 

with Ei = ylf + mf and cr = +1 for a particle and —1 for an antiparticle. In Eq. the 
symmetry factor is denoted by S^y, this will be discussed in the following subsection. 

Next consider the diagrams of Fig. |2|(b) shown schematically in the upper row and in 
more detail in the lower row. To express these diagrams in terms of the scattering operator, 
we define an operator T k such that 

{k, {li}\%c\{Pi}) = (k, {li}\1~\{pi}) — (disconnected parts involving 5{k — Pi)) 

and 

({pi}\T k \k, {/f}) = ({pi}\T\k, {l?}) — (disconnected parts involving 5{k —ft)) . (43) 
We also define 

A k = T-T k . (44) 

The coefficients originating from expanding a cut diagram in this second group must also 
contain a cut with the same restriction that the two external vertices should always be in 
the unshaded region. Hence, the coefficient functions are the zero temperature Cutkosky 
diagrams with the same restriction. We know that the zero temperature Cutkosky diagrams 
represent matrix element of T^T. Hence, the expansion of this group of diagrams must be 

s? ) (*) = i E «— / Il rfr T ^({in\^\{Pi})(k,{Pj}\r\k 7 {t!})\^ 

n>2,a ^{lf} J i=l { Pj } 

+ (disconnected diagrams) 

1 r n-l 

= < E o— / n^<M^}|A&|*,{«f}>U. 

n>2,a J i=l 

+ (disconnected diagrams) . (45) 

The third group of diagrams contain a cut that separates the external vertices as shown 
in Fig. 0(c). Noting that the external frequency k flows out of the unshaded region and the 
cut-line frequencies flow into the shaded region, we can write 

s${k) = i e «— / Ii dr i E(K CT }i^i^fe})^te}|T|{/r})| conn 

n>2,a ^{lf} J i=l { Pj } 

+ (disconnected diagrams) 

. n-l 



1 pn-i 

E *— / I[dr°(k,{i?}\AlA k \k,{in)l 



n>2, CT J i=l 

+ (disconnected diagrams) . (46) 
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Ignoring the disconnected diagrams which will later be shown to cancel, the three con- 
tributions yield 



E* et (£;) = Z^\k) + ZP(k) + E$(h) 

n-l 



= - E -J- I n «flT <*. Of }l (r - iA\r) \k, {zf }> conn 

n>2,a b {lf} J i=l 

= - E O— / II ^ <*> Of {Zf }> conn , (47) 

n>2,,7 ^{If } 3 i=l 

where we used (fc, {Zf }|T fe |A;, {Zf}) conn = (fc, {Zf }|T|fc, {Zf }) conn , Eq. (g|) and the identity 
T = T^S = T\l + iT). 

We can now show that Eq. ([!]) is the lowest order approximation of the expansion (|47|). 
Without self-energy insertions, the matrix element (k,l\A\T\k,l) conn vanishes provided that 
the momentum Z represents a stable particle which cannot decay at T = 0. The lightest 
particle in a theory is certainly stable. Hence, Eq. ([[J) holds at the lowest order in the density 
expansion. 

To check the consistency of Eq. fl4"7p, consider the imaginary part, 



I | . n— 1 

lm£* et (*0 = - ? E «— / II % Of}l Nr k - A{A k ) \k, {inUnn . (48) 

Here we used the unitarity condition ImT = \T^T to get Eq. ( [4~8|) from the second line 
of Eq. (pETD - The cut diagrams leading to expression (|48| ) must contain cuts that separate 
the two external vertices. Also the diagrams corresponding to T^lk must have the external 
momentum k entering the unshaded region, while those corresponding to A^A^ must have 
the external momentum k exiting the unshaded region. In Ref. []I2J it is shown that these cut 
diagrams are exactly the diagrams that gives the spectral density of the 2-point correlation 
function. The retarded self-energy is the negative of a 2-point retarded correlation function 
with a real spectral density (see Appendix [S]). This implies that 

ImS Rct (A:) = -ix s (fc), (49) 



where XsO^) is the spectral density for the self-energy. This is consistent with Eq. ([48]). 



Eq. ( f4"T|) can quickly be used to make contact with the thermal (T > 0) part of the 
thermodynamic grand potential; for a more rigorous discussion see Appendix [D] We need 
to close off the diagrams of Fig. [|by including a propagator for the external line k. The zero 
temperature part will give a contribution that can be included in the scattering matrices, 
while the finite temperature part gives an additional dT integration and the symmetry 
factor must be appropriately adjusted. We note that in Eq. (|47D T^S = J2m=o • 
When the propagator for k is included, 7^ is no longer distinguished from T\ but we must 
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include a factor of {m + 1) 1 in the sum in order to avoid multiple counting. Thus we have 
— % J2m=o (*^) / (m + 1) — i\n(l — iT^ . Thus the thermal part of the grand potential 



can be written 

<>/ . v _L 

n>2,a J i=l 



= -i E «— / n^r({/niin(i-zTt)|{/.r})con„ 

1 n>2,a } J i=l 

1 /• - 

E «— / Il^(R}|ln(l + ^)|{/a)conn, (50) 



where V denotes the volume of the system and in the last step we have used the reality of 
Q. We have also assumed that the physical mass of the particles is used in the propagators 
so that diagrams with a single thermal weighting are superfluous. Equation (]50f) agrees with 
Norton's version [23| of the expression originally given by Dashen, Ma and Bernstein PT . 



A. Symmetry factors 

There are two issues involving the symmetry factor, Sy?}, arising in the equations above. 
One involves the symmetry factor due to the self-interaction. For instance, the two loop 
diagram in a A0 4 theory shown in Fig. ^| has the symmetry factor of S2-100P = 3! since the 
three internal lines are equivalent to each other. Our density expansion inevitably reduces 
the symmetry of this diagram. One needs to show that the expansion in the number of 
automatically generates the symmetry factor associated with the reduced symmetry in such 
cases. The second issue involves identical particles. When there are m identical particles in 
the final state, the m-body phase-space includes a factor of 1/m!. Our expansion must also 
automatically account for this factor. 

To see that the symmetry factors are correctly generated, we label each line in a diagram 
by l a and each vertex by r a and regard them as distinguishable. The symmetry group G 
of a given Feynman diagram then represents the permutations of the internal lines and the 
vertices that does not change the shape of the diagram |24] . The order qc of this group is the 



symmetry factor for the diagram. For example, the symmetry group of the two-loop diagram 
in Fig. [|is the full permutation group of the three identical lines. Hence, S2-I00P = <?g = 3!. 

Consider a cut diagram for the self-energy. The density expansion essentially amounts 
to the sum of all possible ways of mixing the zero-temperature propagators and the thermal 
phase space factors without totally disconnecting the diagram. In such "divided diagrams" 
the lines and the vertices in the original diagram are apportioned into the thermal part 
(which solely consists of T^'s) and the zero-temperature part. 

In general this partitioning will reduce the symmetry group of the original group G to 
a subgroup S. The reduced symmetry group S is then a product group consisting of the 
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FIG. 3. A 2-loop diagram in scalar A</> 4 theory. 



symmetry groups S z and S t of the two parts. Here the subscript 'z' represents the zero- 
temperature part, and the subscript 't' represents the thermal part. In other words the 
symmetry group is now, 

S = S z ® S t . (51) 

The order of this product group is of course qs = q z x q t , where q z and q t are the order of the 
subgroups S z and S t , respectively. We must show that the density expansion automatically 
generates the symmetry factor qs = q z X q t associated with the reduced diagram. 

It is well known in group theory that the order of a subgroup is a factor of the order of 
the full group [25|. That is, qc/Qs — r is an integer. We would like to show that there are 



exactly r different partitions of the original lines and vertices that result in the same divided 
diagram. Then we show that the correct symmetry factor is produced for each individual 
diagram in our expansion. 

To see how it works out, suppose the symmetry group G of the original diagram is of 
order qa, and when divided into zero-temperature and thermal parts, the symmetry group 
reduces to a subgroup S of order q s . Consider the right coset, 

G/S = {Sg u Sg 2 ,---,Sg qG } , (52) 

where each gi G G is a permutation acting on the internal vertices and the lines that does 
not change the shape of the original diagram. Since the identity must be an element of G 
and the Sg^'s are either identical or disjoint, we can rewrite 

G/S = {S 1 ,S 2 ,---,S r } , (53) 



where Si = S, Si's are all disjoint, each Si has qs elements, and r = qc/qs |25| . 

Operating with G Si on the diagram does not change the partition of the lines and 
vertices since Si = S is the symmetry group of the divided diagram. However, the other 
Sj's must contain at least one element that corresponds to a different partition. Then there 
exist hi's such that hi G G and 
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FIG. 4. Density expansion of the 2-loop diagram in Fig. ||. The dot at the end of a line denotes 
that the line corresponds to a thermal particle. 



Si — Shi — {hi, Szhi, s%hi, • • • , s qs hi} . 



(54) 



Here, Si G S, s\ = I and h\ = I. Since Si and Sj are disjoint for i ^ j, we have hi 7^ hj. 
The number of hi 8 is then equal to r = qc/qs- Also each element of Si is related to hi 
by an element of S. This implies that all elements in Si correspond to the same partition 
of the lines and vertices. Hence the total number of inequivalent partitions of the internal 
lines and vertices is exactly the same as the number of h^s. Consequently, when our density 
expansion is made, the coefficient of a divided diagram is r/qa = l/(q z X q t ). Thus the 
symmetry factor associated with the symmetry of the divided diagram is always correctly 
produced by the density expansion. 

As a simple example, consider the 2-loop diagram shown in Fig. |3|. If one of the lines is 
opened out so that it carries a thermal phase space factor as shown in the first of Fig. |], 
the symmetry group is reduced to the permutation of the two remaining lines so that the 
overall factor becomes 1/2!. Originally, qc = 3!, and now q$ = 2!. The ratio qc/qs — 3 of 
course corresponds to the three choices we can make when selecting a line to replace. 

If two lines from the diagram in Fig. are opened out and carry thermal weightings, 
we produce the second diagram in Fig. [§. The symmetry group of the reduced diagram is 
the trivial identity. Hence, the symmetry factor associated with the reduced diagram is 1. 
However, since there are only 3C2 = 3 choices of picking a pair of lines, the overall factor 
for the second diagram in Fig. |4] is only 3/3! = 1/2!. This is as it should be because the 
thermal lines have a reflection symmetry. 

To make the connection between opened-up self-energy diagrams, such as those in Fig. |], 
and the T-matrices label the latter with 2 fc's, m initial momenta Pi and m final momenta 
qi. Eventually, we identify Pi = qi and integrate over the thermal phase-space of all p^s 
to make the contribution to the self-energy. Suppose that all pi and belong to the same 
species. Then the matrix elements 



T(k, p ai , • • • , Pa m \ k, q ai , ■ ■ ■ , q<j m ) , 
where a is one of the m! permutation of {l,...,m}, correspond to the same self-energy 



T(k,px, ...,p m ] k,qi, . 



Am) 



and 



17 



FIG. 5. Cancelling propagators. 



diagram when the p^s are integrated over the thermal phase space. 

Even though there are m! such permutations, the number of distinct T(k, {p}; k, {q}) 
need not equal ml. For some permutations, the above two expressions may not only give the 
same contribution to the self-energy, but actually be the same even before the integration. 
If we reconnect the pt vertices with the vertices, the resulting diagram is a self-energy 
diagram. Hence the number of permutations giving rise to the same T must be equal to 
the order q t of the symmetry group St of the reconnected part of the diagram. Hence, the 
coefficient of this reconnected diagram is given by m\/(q t x q z ) where q z is the symmetry 
factor associated with the T. 

Comparing this with l/(q t x q z ) which we had for the opened-up self-energy diagram, we 
see that a diagram generated from the self-energy is smaller by a factor of m\ than an actual 
scattering amplitude diagram. The generalization to the many species case is immediate. 
Thus if each particle species i has external thermally-weighted lines, the overall factor 
for the opened-up self-energy diagram is = l/flj^i- 



B. Disconnected Parts 



Suppose that, when the internal lines are opened out and given thermal weightings T^, we 
get a disconnected piece containing at least one internal vertex. As long as this disconnected 
piece does not contain the largest time vertex, all cuts for this diagram, including the case 
with no cuts at all, are possible with the rest of the original diagram being kept common. For 
the isolated subdiagram, this situation corresponds to having all the frequency denominators 
in Eq. (|28|) and summing over all time orderings. Consequently, the net contribution of such 
a disconnected part vanishes. In the real-time method this cancellation is referred to as the 
vanishing of the sum of all circlings ||. 

For an elementary example, consider the sum of all cut and uncut propagators depicted 
in Fig. [|. The sum is 

(-ifG c (k) + (i) 2 G* c (k) + (-i){i)A^(k) + (z)(-i)A+(k) = , (55) 

where all the propagators are zero-temperature ones. From the structure of the propagators 
(j35p - (p8|), it is clear that the sum is always zero. 
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FIG. 6. Schematic depiction of the possible ways to cut a self-energy insertion diagram. 



C. Self-Energy Insertions 



For simplicity we only consider scalar particles without a chemical potential. The argu- 
ment presented here can be generalized immediately to other cases. Self-energy insertions 
are potentially hazardous when the propagators have poles at all four k° = ±Ek ± is. This 
is the case for us due to the addition of the thermal phase space factor to the usual zero- 
temperature propagators (c.f. Eqs. (|35|) - fl38f )). The product of any two propagators with 
the same argument then contains pinching poles. 

To show that in fact the pinching pole contributions all cancel, consider the diagrams 
shown in Fig. They differ only in the manner of cutting. Representing the blobs in the 
first and third diagrams by C(k) and L(k), respectively, the sum of the four diagrams can 
be written 

F{k) = G B (k)C(k)G B (k) + A B (k)C*(k)A+(k) 

- A B (k)L(k)G B (k) - G B (k)L(-k)A+(k) . (56) 

Extracting the coefficient of the pinching poles 2(k 2 — m 2 + ie)~ l {k 2 — m 2 — ie)^ 1 , we find 

^pinchW = n B (k°)N B (k°) (C{k) + C*(k)) 

- n B (fco)(l/2 + n B {k ))L{k) - N B {k )(l/2 + n B (k ))L(-k) . (57) 

If .Fpinch(fc) were non-zero, a self-energy insertion would cause an uncontrollable divergence. 
Fortunately, J r P m C h(k) does vanish due to the following properties of C(k) and L(k) p|,|7|,|8|JT6[| 

C(k) + C*(k) = L(k) + L(—k) (58) 

and 

L(-k) = e- k0f3 L(k) . (59) 

The first identity (|58|) is the finite temperature version of the optical theorem. The second 
identity can be easily obtained by using A B (k) = e k0/3 A B (k). Therefore, pinching poles do 
not occur in self-energies insertions. 

The absence of pinching poles has two important consequences. First, self-energy in- 
sertions do not cause uncontrollable divergences due to the cancellation between cut and 
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uncut self-energies. Second, the propagators connected to a self- energy do not produce ther- 
mal phase-space factors. Without this second point it would not be possible to identify the 
coefficients of the expansion in the factors with the scattering amplitudes. 



D. Polarization Factors and the Overall Sign 

To convert a Feynman diagram to a scattering amplitude, one needs the polarization 
factors for the external lines. For an external gauge boson line, a polarization vector e M (p, s) 
orthogonal to the incoming momentum p is needed. This is provided by the substitution 

9nv^~ e »(P, s ) e t(P, s ) ( 60 ) 

s=l,2 

inside a Feynman diagram; this is valid due to the Ward identity. In the Feynman gauge, 
the gauge boson phase space factor is proportional to the metric g^ u which then provides 
appropriate e^'s. 

For fermions, each incoming fermion (anti-fermion) line requires a factor of the Dirac 
spinor u s (p) (v s (p)), and each outgoing fermion (anti-fermion) line requires a factor of u s (p) 
(v s (p)). For the external line corresponding to the self-energy momentum, one can use the 
Dirac spinor identity 

1 = ^ E Mk)u s (k) - v s {k)v s {k)) , (61) 
where the spinors are normalized according to Peskin and Schroeder [2(|]. This yields 



E F (Jfe) = l£ F (fc)l 

= -t—~ Y\u s (k)\u s (k)T, F (k)u s (k)]u s (k) + v s (k)\v s (k)'Ep(k)v s (k)\v s (k)} . (62) 
Arrip g 1 L J L J > 

The factor u s (k)Y,F(k)u s (k)/(2mF) then has a multiple scattering expansion in terms of 
spin up-up and down-down scattering amplitudes involving the particles, while the factor 
v s (k)T>F(k)v s (k)/(2mF) has a multiple scattering expansion in terms of spin up-up and 
down-down scattering amplitudes involving the anti-particles. The mixed term vanishes 
since u(k)j°v(k) = v(k)^°u(k) = u(k)k ,I 'j l j,v(k) = v(fc)fc /i 7 At w(fc) = 0, and the self-energy, as 
well as the spectral density, must have the structure 

Ep(ib) = A(T, fc)7% - B(T, k)m F + C(T, k)^ (63) 

at finite temperature. 

As an example, consider the nucleon self-energy in a thermal pion medium. Considering 
only the strong interaction, we can regard both the pions and the nucleons as stable. The 
up-up component of the lowest order nucleon self-energy can then be expressed as 
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FIG. 7. The one-loop electron self-energy. The curly line represents the photon. 
u + (/c)St(A;)m + (A;) 



S++(A;) 



2,771 TV 
1 „ f d% 



E / ( 2n )32E T N^Nn a (k + l n ^k + Q 



2m N ^J (2tt)32^ 

d 3 ^ /r, x ( Vs 



where T^ a ^ Nna (k + l n —> k + l n ) is the sum of all scattering amplitudes including the cross 
terms. In the last line we used T = 8iry/sf. This expression differs by a factor of y/s/m^ 
from that given by Eletskii and Ioffe ||26|| . At low temperature (T < ffljv) and in the nucleon 



rest frame, their expression may be justified because the pion thermal energy is insignificant 
compared to the nucleon rest mass. 

For internal fermion lines which are opened up and given thermal weightings we use the 
relation 

(F 7m + m F ) 2tt 5(k 2 - E 2 k ) = (A£ 7a , + m F ) ^ 5(k - E k ) - (A£ 7 „ - m F ) ^ 5(k + E k ) 

2n 2n 
= Y, u s(k+)u s (k + ) — - 5(k - E k ) — E v s (k_)v s (k_) — — 5(k + E k ) , (65) 

s s 

where k± = (E k ,±k). Hence for both spin-1 gauge bosons and spin-| fermions the polar- 
ization factors are all correctly accounted for. 

When fermions are involved, our expansion must generate the correct overall sign of a 
diagram. The presence of a fermion loop in a diagram carries an additional overall factor of 
(— 1). When the fermion loop is broken this sign is carried by the thermal phase space factor, 
Tp rather than the scattering amplitude. The vv part of Eq. QB5| ) carries an additional (—1) 
due to the exchange of the fermion legs. When a fermion line which is not a part of a fermion 
loop carries a thermal phase space factor, the uu term corresponds to a crossed diagram and 
the required factor of (—1) is provided by Tp. Hence, the overall signs are also correctly 
accounted for in our expansion. 



IV. ELECTRON SELF-ENERGY 

As an example of applying the method developed here, consider the electron self-energy 
at temperatures T <C m e obtained from the diagram of Fig. [?]. At these low temperatures, 



21 



we can neglect the electron thermal correction because it is suppressed by a factor of e~ l3me . 
The self-energy written down for the nucleon ( S4|) is valid for the electron if the pion is 
changed to a photon with the summation referring to the photon polarization: 

d?k 



2m, 



E 



e s =l,2 



(27r) 3 2|k| 



n B (\k\)T++ (p + k^p + k) . 



(66) 



After the photon polarization summation, the spin up-up component of the tree-level Comp- 
ton scattering amplitude is given by (e.g. Ref. |20| ) 



r £t ei {v + k -> p + k) = e 2 u+(p) 
In the electron rest frame, this reduces to 
T++ n (p + k^p + k) = 
The real part in the rest frame is then, 



2p-k + ie 



+ 2j) | - 7 ^ 7/ , + 2/ 
—2p-k + ie 



u+(p) 



(67) 



-4e 2 



1 + 



2m? e ei 



Am 2 p \k\ 2 + e 2 



(68) 



2e 2 

KeHt + (m P ) = 



d 3 k 



m e J (2vr) 3 2|k| 
e 2 T 2 vraT 2 



12m P 



3m» 



(69) 



The result, of course, coincides with a previous calculation ||27|| . The imaginary part is 

d 3 k n fl (|k|) 



Im (m e 



lim 4e 2 m«e 



(27r) 3 2|k| (4m 2 |k| 2 + e 2 ) 



47T 



(70) 



This agrees with the leading term given by Henning et al. EH . 



V. CONCLUSION 



In this paper, the multiple scattering expansion of the thermal correction to the retarded 
self-energy is presented starting from the imaginary-time formalism. The leading order term 
of this expansion corresponds to the often used rule that the thermal correction to the self- 
energy is the thermal phase-space times the scattering amplitude. Although the formal 
expansion can be always made, we have argued that the expansion is useful only if the 
minimum of \m a ± fi a \ is large compared to the temperature. We have also demonstrated 
the connection between the self-energy and the thermal part of the grand potential. 

The result presented here may be used in two ways. One is to use existing calculations 
of scattering amplitudes to calculate the self-energy. In this way, the considerable effort 
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usually needed to evaluate thermal correlation functions can be much reduced. The other 
way is simply to use experimental scattering amplitudes to calculate the self-energy. For 
interactions involving large coupling constants, this may be the only reliable way to calculate 
the thermal correction to the self-energy. 

The reliable calculation of medium effects is important in analyzing data from heavy- 
ion collision experiments. For instance, how the p-meson behaves in-medium can greatly 
influence the dilepton spectrum in heavy ion collisions [p],p9| . The in-medium effect becomes 
even more important in future RHIC (Relativistic Heavy Ion Collider) experiments. The 
main goal of RHIC is to find the quark-gluon plasma. Hence, it is crucial to understand the 
hadronic part of the in-medium finite-temperature effect so as to separate this signal from 
that of the long-sought quark-gluon plasma. Even though a full theory of hadron scattering 
is lacking, there is a considerable amount of data on scattering cross-sections accumulated 
in the past decades. Utilizing those data we may at least phenomenologically separate a 
truly hadronic effect from the effect of a new state of matter. 
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APPENDIX A: EUCLIDEAN PROPAGATORS 



1. Bosonic Propagator With Chemical Potential 



When the chemical potential is non-zero, we must at least deal with a complex field. 
The effective Hamiltonian is 

k = H-^2^aQa. (Al) 

a 

For notational convenience we suppress spatial indices in this section. 
The spectral density for the propagator is defined to be |11 

m,n 

= ^5(u + K m - K n ) (e~ 0Km - e-^») mn </>L 
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(l - e~^) ^S(u + K m - K n ) 



The Euclidean propagator is given by 

G b (t) = (T0(r)0t( O )) 

= 6(t)Ti (e-^ k e Tk 4>e- Tk ^) + 9(-r)Tr (f^ <fie rk $e~ T 

= \% (e-^[l + n B {uj)}p B {uj)6{T) + e-" T n B {uj)p B {uj)6{-T 
= J % [1 + n B (u)\ (e-^p+iu) 9{r) + p B (co) 6(-r 

where we used 



tib(-uj) = -[1 +n B (u)} 



and defined 



Pb{u) = p B (uj) ; Pb (uj) = -p B {-uJ) . 
To find p B (w), use (|TT| 

G B (iu B ,k) = — 2 ■ — ^ = / : 

E-f. — [iv B + p) 1 J 2n uj — iv B 



with Ek = Vk 2 + m 2 , and take the discontinuity across the real axis to get 
p^(cj,k) = sign(u; + p) 2tt5 ((u; + /i) 2 - E 2 ^ 



and 



p B (a>,k) = -p£(-a;,k) = sign(w - p)2n5 ({uj - p) 2 - El) . 



2. Fermionic Propagators with Chemical Potential 

The spectral density for the fermion propagator is defined to be 



£ 2<k5(u + K m - K n ) (e-^™ + e"^") z/wA, 

m,n 

(l + e^) £ 27r5(o; + K m - tf n ) VwVWne~^ 



where the brace denotes an anticommutator and we used the fact that matrix elements 
as ip mn are c-numbers. 

The Euclidean propagator is given by 
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G F (r) = (TV>(t$(0)> 

= 6{t)Ti {e- pk e Tk il)e- Tk f) - 0(-r)Tr (e^ ^e Tk ^e~ Tk 

% (e— [1 - n F ( W )] 4(o;)fi(r) - n P ( W ) p+ M(9(-r)) 
^ [1 - n F (u,)] (e— p+H^(r) + e -p F (^(-r)) , 
where the change of variable uo — > — c<j was made for the second term and we used 

n F (-a;) = ^ + 1 = 1 - n F M . 

Also, we defined 

PfM = PfM ; pFM = -pf(-w) ■ 

Hence, for both bosons and fermions we have 

G c (r)= f^N c (u)(e— p c >)0(r) + p c ~(u;)0(-t)) , 



with 



and 



N c {u) = 1 + (-l) C n c (a;) = 0(u) + (-l)<sigii(u;)n f (M) , 

where 

(-1) B = 1 and (-1) F = -1 . 

To find pp(u), start from the Euclidean propagator 

z (V F - z» 7° - p-7 + m _ f <ko p F (o;,k) 
<^F^F,PJ - 7 . u , n, - / 7, : > 

with Ep = y/p' 2 + m' 2 , and take the discontinuity across the real axis to get 
f4(u, p) = ((u + /i)7° - P-7 + m) sign(^ + p) 2-k5((u + p,) 2 - E 2 ) . 



APPENDIX B: REAL TIME PROPAGATORS 



For both bosons and fermions the real time version of Eq. ( [A 1 3| ) is 



2vr 
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The Fourier transform yields the momentum space propagator 



p [»(«) + (-l)<sign(u,H(M)] -AW + f"' , (B2) 



using Eq. ( A15 ). The propagator has two terms. The first corresponds to the zero temper- 
ature case 



Qg(*)=/gW P ? U) , ,) , (B3) 

which yields the standard Minkowski propagators given in the text. The second term cor- 
responds to the finite temperature phase space factor 



r C (*) = (-D £ / ^sign(uOn c (M) ( + — £ 



2vr ° w ' Vl ye + i(a> - A; ) e + i(u; + A; ) / 
= (-l)^ign(k )n c (\k°\)p+(k) , (B4) 

where we used p^{uj) = —p^(—u). 

Cut propagators can also be decomposed into zero and non-zero temperature parts. The 
cut propagator for a particle is 

A+(k) = N c (k°)p+(k) = 6(k°)p+(k) + (-l)«sign(fc°)n c (|A; |) p+(k) . (B5) 

Letting the momentum and frequency follow the charge, we get the cut propagator for an 
anti-particle N^(k°)p^(k) — > N^(—k°)p^(—k) which is 

A f -(fc) =N ( (-k°)p<(-k) = (-l)<n c (k°)p+(k) 

= -0(-k°)p+(k) + (-l)i sign(fc°) n f (|A; |) , (B6) 

using 

n c (A; ) = -(-l)<0(-fc°) + sign(k )n c (\k°\) . (B7) 
We also note that Af(fc) = (-l) c e^ fco A^(A;). 

APPENDIX C: DERIVATIVE COUPLINGS 

For the time derivative of the propagator in Eq. (|A13|) we have 
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|:G c (r,k) = J ^N ( (lu) (-cop+(u,k)e-^9(r) + up^u, k) e" T 9(-r) 
+ J f^c(«0 (pf(o;,k)-pf ( W) k))*(r) 
^iV c H (-o;p+(a;,k) e -^(r) + c^(u,, k) e^0(-r) 

+ /^pf(a;,k)5(T). (CI) 

Thus the spectral density for the time derivative of the propagator is —up^(uj, k) and there 
is also a 5(t) term. The latter vanishes for bosons since / cL;pg(u;,k) = 0. In general, we 
can say 

d r T G c (r,k) = G c (r,k) + X>«(T)F,(k) , (C2) 

z=o 

where the label £ includes both ( and the number of derivatives r so that has the 
spectral density (— uj) r p^. In the second term of (|C2| ) 5^(r) = <9[<5(r) and the sum does 
not contribute for r = 0. Then, after all the 5-functions in Eq. ( |C2| ) are integrated over, we 
have 

C { p +1 ({qi,iui}) 

rP v ( N \ r d z kr 



i=o \ i=o / " Ler ae r 

+e/ n^exp / n 7o4 n 

+e / n*»«p M / n j^A 2 aK,iv r ,iu q }) n c^w-^) 

g,r - 70 ^ i^ r>g / J Ler \ Z7T ) a&V" 

+ • • • , (C3) 

where in V a pair of vertices refer to the same time, in V" two pairs of vertices refer to 
the same times, and so on. The bars over the V\ indicates that some of them are now 
combinations of the original zVs. 

If a total of n derivatives is contained in the expression for the diagram T, then performing 
the time integrals we have 

n 

tf#+i({q*,«4}) = E E cSH({cu,m}) (C4) 

r=or^cr r 

where 

=/n|ln(/ ^N CaM Pt(k a )) A r ({Kh{m}) n ( a j- 



Lgr V 7 QgT r \ / intervals 

V-r>j>l 



(C5) 
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We know that in the (3 — > oo limit, the frequency denominators we get from Eq. ( |C5|) must 
add up to make the Wick rotated d^G^(t). That is, the coefficients A r must be such that if 
we change iV — > 6, ivi — > qf and add — it to AJ in Eq. ( 05 ), we get the zero temperature 
iV + 1-point function: 

csimum}) - D { p +1 ({ qi }) = i v r n n g c«(^) , (C6) 

as in Eq. ([Z7|). Arguments similar to those in Subsec. [11 B| can be applied to cut diagrams so 
that summing the time-ordered diagrams on the unshaded side leads to the usual Feynman 
rules, while for the shaded side complex conjugate Feynman rules apply. Thus the multiple 
scattering expansion of the self-energy in Sec. 3 is also applicable when derivative couplings 
are present. 



APPENDIX D: THE THERMODYNAMIC POTENTIAL 

Here we consider the multiple scattering expansion of the thermodynamic grand potential 
defined by 

^L = ^J_l nZ (Dl) 

V [3V K J 

where Z is the partition function and V is the volume. In the imaginary-time formalism, 
the thermodynamic potential is the sum of all connected "vacuum" graphs. Analytic con- 
tinuation of such a result may at first appear to be a poorly defined concept since there are 
no external frequencies to start with. Nevertheless, it is possible to consider the vacuum 
graphs as the zero frequency limit of the Appoint functions []30| and, further, it is sufficient 
to consider retarded correlation functions due to the reality of the thermodynamic potential. 

To calculate Q, we regard the (p in Eq. ( PUD as "external" interaction vertices which will 
contain several fields. The retarded functions then consist entirely of all possible interaction 
vertices and there are zero external frequencies and momenta entering or leaving the diagram. 
Eq. (P3"|) can then be written 



x 



n a j+^ n A /-* e • (D2) 



intervals intervals 
V>j>v+1 v>j'>l 

Here the fixed vertex v is arbitrarily chosen for each diagram. Instead of the identity 
we use 

1 1 



AJ + ie A] - ie 



- 2<Ki5(A°) , (D3) 
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FIG. 8. Diagrams for the thermodynamic potential. The dashed line indicates a cut. 



to get 



*fti(o) = e / n S n (7 ^ t w/«w) ^ 

r CT cr J Ler \ zn ) «er V z/l / 



— le 



n 

intervals 
V>j'>l 

+h) £ 2tt5(aj) n ( A fc- 



y>fc>i 



+ H) 2 £ 27r < y(AJ)27rtf(AJ) II (A\ 



ye 



V>«>1 

Ma,* 



+ ... + (-« 



.V— u 



n 2tt5(aj) n ( a ^- 



intervals 
V>j>tH-l 



v>k>l 



-i 



(D4) 



In this way the denominators contain only — ie and consequently, after summing over all 
time orderings a, the result can be expressed entirely in terms of G^(k) and A^(k), i.e. the 
complex conjugate propagator G^(k) does not occur. The price paid is the appearance of 
multiple cuts since each of the ^-functions in Eq. ( p4| ) represents a cut. So the result of 
the time ordering summation will yield uncut diagrams plus those with a sequence of cuts 
illustrated in Fig. |j. The blobs in these diagrams involve the propagators G^(k), while the 
cut lines require according to the direction of the momentum, as before. 

Following the same procedure as before, we can now expand the diagrams in the number 
of thermal phase space factors T^. The contribution of the first term in Eq. ( |D4T ) to the 
thermal part of Q is 
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q(!) i r n 

V n>2,a b m J i=l 

where we have excluded diagrams with a single thermal weighting since we assume that the 
physical masses of the particles are used in the propagators. The contribution of diagrams 
with (m — 1) cuts is 



(m) 



1 r n 

- E c— / li dr ° ({in\%H , rr- i \{inUn . <m 

Here T v contains the arbitarily-chosen fixed vertex, v. If, however, we remove this restriction 
and allow v to lie in any of the T matrices, while compensating for the multiple counting, 
we obtain 



. ^ l r£r _ l 



E O— / Il^-({^}l(-^) m |{ / r})conn • (D7) 



Summing over the number of cuts, we then have 

Q oo ry(m) , n -l 

1 7 : =E^ = ^E^/n^ <{if}| ln(l + iT)\{in) CO n» ■ (D8) 

V m=l V n>2,a b m J i=l 

This is the expression given in Eq. ([5(]) and, as we remarked, it agrees with the result 
Norton [2^] obtained using a different approach. This in turn is equivalent to the expression 
obtained long ago by Dashen, Ma and Bernstein |H from a non-relativistic analysis. 
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